Applicability of leaf area index products for boreal regions of Sweden

Leaf area index (LAI) of boreal ecosystems were estimated with optical instruments at the Laxemar and the Forsmark investigation areas in Sweden. The aim was to study the relationship between LAI and the normalized difference vegetation index (NDVI) from Landsat-5 and SPOT and evaluate the applicability of the MODIS (Moderate Resolution Imaging Spectroradiometer) LAI product for small boreal regions. Relationships between ground-estimated LAI and NDVI were significant for coniferous, deciduous and mixed forest sites in Laxemar. For Forsmark, effective LAI was correlated to NDVI for all sites. LAI estimated from NDVI was also used for evaluating accuracy of the MODIS LAI product. The comparison showed no correlation between MODIS LAI and NDVI-based LAI in Forsmark whereas there was in Laxemar. MODIS LAI was on average 2.28 higher than NDVI-based LAI and it also showed larger scatter. Scale issues were the main explanation to high MODIS LAI, since the heterogeneous landscapes with open areas (given a value of zero in the NDVI estimates) was seen as forest in the large pixels of the MODIS LAI product. Therefore, we do not recommend using the MODIS LAI product in small boreal regional landscapes, such as the Forsmark and Laxemar investigation areas

Comparing parametric and non-parametric approaches for estimating trends in multi-year NDVI

The aim of this study is to systematically compare parametric and non-parametric techniques for analyzing trends in annual NDVI derived from NOAA AVHRR sensor in order to examine how trend type and departure from normality assumptions affect the accuracy of detecting long-term change. To generate annual data, the mean NDVI of a four-month long ‘green’ season was computed for fifteen sites (located in Africa, Spain, Italy, Sweden, and Iraq) from the GIMMS product for the periods 1982-2006. Trends in these time series were then estimated by Ordinary Least-Squares (OLS) regression (parametric) and the combined Mann-Kendall test with Theil-Sen slope estimator (non-parametric), and compared using slope value and statistical significance measures. We also estimated optimal polynomial model for the annual NDVI, by using Akaike Information Criterion (AIC), to determine the trend type at each site. Results indicate that slopes and their statistical significances obtained from the two approaches at sites with low degree polynomials (mostly linear) and steep monotonic (gradually increasing or decreasing) trends compare favourably with one another. At sites with weak linear slopes, the two approaches had similar results as well. Exceptions include sites with abrupt step-like changes resulting in departures from linearity and consequently high degree polynomials where the least-squares method outperformed the Mann-Kendall Theil-Sen method. In sum, we conclude that OLS is superior for detecting NDVI trends using annual data though further investigation using other techniques is recommended

The supply and demand of net primary production in the Sahel

Net primary production (NPP) is the principal source of energy for ecosystems and, by extension, human populations that depend on them. The relationship between the supply and demand of NPP is important for the assessment of socio-ecological vulnerability. We present an analysis of the supply and demand of NPP in the Sahel using NPP estimates from the MODIS sensor and agri-environmental data from FAOSTAT. This synergistic approach allows for a spatially explicit estimation of human impact on ecosystems. We estimated the annual amount of NPP required to derive food, fuel and feed between 2000 and 2010 for 22 countries in sub-Saharan Africa. When comparing annual estimates of supply and demand of NPP, we found that demand increased from 0.44 PgC to 1.13 PgC, representing 19% and 41%, respectively, of available supply due to a 31% increase in the human population between 2000 and 2010. The demand for NPP has been increasing at an annual rate of 2.2% but NPP supply was near-constant with an inter-annual variability of approximately 1.7%. Overall, there were statistically significant (p < 0.05) increases in the NPP of cropland (+6.0%), woodland (+6.1%) and grassland/savanna (+9.4%), and a decrease in the NPP of forests (−0.7%). On the demand side, the largest increase was for food (20.4%) followed by feed (16.7%) and fuel (5.5%). The supply-demand balance of NPP is a potentially important tool from the standpoint of sustainable development, and as an indicator of stresses on the environment stemming from increased consumption of biomass

Classification of grassland successional stages using airborne hyperspectral imagery

Plant communities differ in their species composition, and, thus, also in their functional trait composition, at different stages in the succession from arable fields to grazed grassland. We examine whether aerial hyperspectral (414–2501 nm) remote sensing can be used to discriminate between grazed vegetation belonging to different grassland successional stages. Vascular plant species were recorded in 104.1 m2 plots on the island of Öland (Sweden) and the functional properties of the plant species recorded in the plots were characterized in terms of the ground-cover of grasses, specific leaf area and Ellenberg indicator values. Plots were assigned to three different grassland age-classes, representing 5–15, 16–50 and >50 years of grazing management. Partial least squares discriminant analysis models were used to compare classifications based on aerial hyperspectral data with the age-class classification. The remote sensing data successfully classified the plots into age-classes: the overall classification accuracy was higher for a model based on a pre-selected set of wavebands (85%, Kappa statistic value = 0.77) than one using the full set of wavebands (77%, Kappa statistic value = 0.65). Our results show that nutrient availability and grass cover differences between grassland age-classes are detectable by spectral imaging. These techniques may potentially be used for mapping the spatial distribution of grassland habitats at different successional stages

Radiometric Correction of Multispectral UAS Images: Evaluating the Accuracy of the Parrot Sequoia Camera and Sunshine Sensor

Unmanned aerial systems (UAS) carrying commercially sold multispectral sensors equipped with a sunshine sensor, such as Parrot Sequoia, enable mapping of vegetation at high spatial resolution with a large degree of flexibility in planning data collection. It is, however, a challenge to perform radiometric correction of the images to create reflectance maps (orthomosaics with surface reflectance) and to compute vegetation indices with sufficient accuracy to enable comparisons between data collected at different times and locations. Studies have compared different radiometric correction methods applied to the Sequoia camera, but there is no consensus about a standard method that provides consistent results for all spectral bands and for different flight conditions. In this study, we perform experiments to assess the accuracy of the Parrot Sequoia camera and sunshine sensor to get an indication if the quality of the data collected is sufficient to create accurate reflectance maps. In addition, we study if there is an influence of the atmosphere on the images and suggest a workflow to collect and process images to create a reflectance map. The main findings are that the sensitivity of the camera is influenced by camera temperature and that the atmosphere influences the images. Hence, we suggest letting the camera warm up before image collection and capturing images of reflectance calibration panels at an elevation close to the maximum flying height to compensate for influence from the atmosphere. The results also show that there is a strong influence of the orientation of the sunshine sensor. This introduces noise and limits the use of the raw sunshine sensor data to compensate for differences in light conditions. To handle this noise, we fit smoothing functions to the sunshine sensor data before we perform irradiance normalization of the images. The developed workflow is evaluated against data from a handheld spectroradiometer, giving the highest correlation (R-2 = 0.99) for the normalized difference vegetation index (NDVI). For the individual wavelength bands, R-2 was 0.80-0.97 for the red-edge, near-infrared, and red bands

Development of a method for monitoring of insect induced forest defoliation - limitation of MODIS data in Fennoscandian forest landscapes

We investigated if coarse-resolution satellite data from the MODIS sensor can be used for regional monitoring of insect disturbances in Fennoscandia. A damage detection method based on z-scores of seasonal maximums of the 2-band Enhanced Vegetation Index (EVI2) was developed. Time-series smoothing was applied and Receiver Operating Characteristics graphs were used for optimisation. The method was developed in fragmented and heavily managed forests in eastern Finland dominated by Scots pine (Pinus sylvestris L.) (pinaceae) and with defoliation of European pine sawfly (Neodiprion sertifer Geoffr.) (Hymenoptera: Diprionidae) and common pine sawfly (Diprion pini L.) (Hymenoptera: Diprionidae). The method was also applied to subalpine mountain birch (Betula pubescens ssp. Czerepanovii N. I. Orlova) forests in northern Sweden, infested by autumnal moth (Epirrita autumnata Borkhausen) and winter moth (Operophtera brumata L.). In Finland, detection accuracies were fairly low with 50% of the damaged stands detected, and a misclassification of healthy stands of 22%. In areas with long outbreak histories the method resulted in extensive misclassification. In northern Sweden accuracies were higher, with 75% of the damage detected and a misclassification of healthy samples of 19%. Our results indicate that MODIS data may fail to detect damage in fragmented forests, particularly when the damage history is long. Therefore, regional studies based on these data may underestimate defoliation. However, the method yielded accurate results in homogeneous forest ecosystems and when long-enough periods without damage could be identified. Furthermore, the method is likely to be useful for insect disturbance detection using future medium-resolution data, e. g. from Sentinel-2.Peer reviewe

Upscaling Northern Peatland CO2 Fluxes Using Satellite Remote Sensing Data

Peatlands play an important role in the global carbon cycle as they contain a large soil carbon stock. However, current climate change could potentially shift peatlands from being carbon sinks to carbon sources. Remote sensing methods provide an opportunity to monitor carbon dioxide (CO2) exchange in peatland ecosystems at large scales under these changing conditions. In this study, we developed empirical models of the CO2 balance (net ecosystem exchange, NEE), gross primary production (GPP), and ecosystem respiration (ER) that could be used for upscaling CO2 fluxes with remotely sensed data. Two to three years of eddy covariance (EC) data from five peatlands in Sweden and Finland were compared to modelled NEE, GPP and ER based on vegetation indices from 10 m resolution Sentinel-2 MSI and land surface temperature from 1 km resolution MODIS data. To ensure a precise match between the EC data and the Sentinel-2 observations, a footprint model was applied to derive footprint-weighted daily means of the vegetation indices. Average model parameters for all sites were acquired with a leave-one-out-cross-validation procedure. Both the GPP and the ER models gave high agreement with the EC-derived fluxes (R-2 = 0.70 and 0.56, NRMSE = 14% and 15%, respectively). The performance of the NEE model was weaker (average R-2 = 0.36 and NRMSE = 13%). Our findings demonstrate that using optical and thermal satellite sensor data is a feasible method for upscaling the GPP and ER of northern boreal peatlands, although further studies are needed to investigate the sources of the unexplained spatial and temporal variation of the CO2 fluxes.Peer reviewe

Cross-sectoral impacts of the 2018–2019 Central European drought and climate resilience in the German part of the Elbe River basin

The 2018–2019 Central European drought was probably the most extreme in Germany since the early sixteenth century. We assess the multiple consequences of the drought for natural systems, the economy and human health in the German part of the Elbe River basin, an area of 97,175 km2 including the cities of Berlin and Hamburg and contributing about 18% to the German GDP. We employ meteorological, hydrological and socio-economic data to build a comprehensive picture of the drought severity, its multiple effects and cross-sectoral consequences in the basin. Time series of different drought indices illustrate the severity of the 2018–2019 drought and how it progressed from meteorological water deficits via soil water depletion towards low groundwater levels and river runoff, and losses in vegetation productivity. The event resulted in severe production losses in agriculture (minus 20–40% for staple crops) and forestry (especially through forced logging of damaged wood: 25.1 million tons in 2018–2020 compared to only 3.4 million tons in 2015–2017), while other economic sectors remained largely unaffected. However, there is no guarantee that this socio-economic stability will be sustained in future drought events; this is discussed in the light of 2022, another dry year holding the potential for a compound crisis. Given the increased probability for more intense and long-lasting droughts in most parts of Europe, this example of actual cross-sectoral drought impacts will be relevant for drought awareness and preparation planning in other regions

Cross-sectoral impacts of the 2018–2019 Central European drought and climate resilience in the German part of the Elbe River basin

The 2018–2019 Central European drought was probably the most extreme in Germany since the early sixteenth century. We assess the multiple consequences of the drought for natural systems, the economy and human health in the German part of the Elbe River basin, an area of 97,175 km2 including the cities of Berlin and Hamburg and contributing about 18% to the German GDP. We employ meteorological, hydrological and socio-economic data to build a comprehensive picture of the drought severity, its multiple effects and cross-sectoral consequences in the basin. Time series of different drought indices illustrate the severity of the 2018–2019 drought and how it progressed from meteorological water deficits via soil water depletion towards low groundwater levels and river runoff, and losses in vegetation productivity. The event resulted in severe production losses in agriculture (minus 20–40% for staple crops) and forestry (especially through forced logging of damaged wood: 25.1 million tons in 2018–2020 compared to only 3.4 million tons in 2015–2017), while other economic sectors remained largely unaffected. However, there is no guarantee that this socio-economic stability will be sustained in future drought events; this is discussed in the light of 2022, another dry year holding the potential for a compound crisis. Given the increased probability for more intense and long-lasting droughts in most parts of Europe, this example of actual cross-sectoral drought impacts will be relevant for drought awareness and preparation planning in other regions

Tree canopy extent and height change in Europe, 2001-2021, quantified using Landsat data archive

European forests are among the most extensively studied ecosystems in the world, yet there are still debates about their recent dynamics. We modeled the changes in tree canopy height across Europe from 2001 to 2021 using the multidecadal spectral data from the Landsat archive and calibration data from Airborne Laser Scanning (ALS) and spaceborne Global Ecosystem Dynamics Investigation (GEDI) lidars. Annual tree canopy height was modeled using regression tree ensembles and integrated with annual tree canopy removal maps to produce harmonized tree height map time series. From these time series, we derived annual tree canopy extent maps using a >= 5 m tree height threshold. The root-mean-square error (RMSE) for both ALS-calibrated and GEDI-calibrated tree canopy height maps was = 94% for the tree canopy extent maps and >= 80% for the annual tree canopy removal maps. Analyzing the map time series, we found that the European tree canopy extent area increased by nearly 1% overall during the past two decades, with the largest increase observed in Eastern Europe, Southern Europe, and the British Isles. However, after the year 2016, the tree canopy extent in Europe declined. Some regions reduced their tree canopy extent between 2001 and 2021, with the highest reduction observed in Fennoscandia (3.5% net decrease). The continental extent of tall tree canopy forests (>= 15 m height) decreased by 3% from 2001 to 2021. The recent decline in tree canopy extent agrees with the FAO statistics on timber harvesting intensification and with the increasing extent and severity of natural disturbances. The observed decreasing tree canopy height indicates a reduction in forest carbon storage capacity in Europe